Fermionic out-of-plane structure of polarization singularities 
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A new classification of circular polarization C points in three-dimensional polarization ellipse fields is proposed. 
The classification type depends on the out-of-plane variation of the polarization ellipse axis, in particular, 
whether the ellipse axes are in the plane of circular polarization one or three times. A minimal set of parameters 
for this classification are derived, and discussed in the context of the familiar in-plane C point classification into 
lemon, star, and monstar types. This new geometric classification is related to the Mobius index of polarization 
singularities recently introduced by Freund. © 2011 Optical Society of America 
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A goal of modern optics is a complete understanding 
of the fine interference structure of classical monochro- 
matic fields. At its most general, a nonparaxial, complex 
vector light field in three dimensions has elliptic polar- 
ization at each point, in part characterized by the ellipse 
major axis, whose direction in three dimensions varies 
with position. The ellipse axis is a director (represented 
by ± a vector), and this director direction is undefined 
and topologically singular when the polarization is circu- 
lar, which occurs along space curves in three dimensions 
called 'C lines' [I]. These C lines are topological defects 
in the smooth complex vector fields describing the light 
field, and on a loop around the C line, the director un- 
dergoes a rotation through n, fermionically exchanging 
its two sides (i.e. switching sign after a rotation by 2ir). 
These polarization singularities have been the subject of 
much study (reviewed in Refs. HJG2), an d 2D properties 
have been measured [4j[5] , although measurement of 3D 
fields is still a challenge at optical frequencies. 

Topological defects have topological and geometric 
properties, and in Ref. |TJ Nye and Hajnal identified 
several of these for points on C lines looking at the 
nearby ellipses projected into the plane of the circular 
polarization like their 2D counterparts [6j: e.g. the ±1/2 
topological index of the director's rotation about the C 
point, and the 'L classification' [7] according to which 
the projected axis becomes purely tangential to a circle 
around the C point once or three times, denoted LI or 
L3. There are three types of two-dimensional C point, 
'lemon' (+1/2 LI), 'star' (—1/2 L3) and the transitional 
'monstar' (+1/2 L3), and these are determined by two 
parameters describing the nonuniform rate of change of 
the director around the C point [5]. 

However, such properties are purely 2-dimensional, 
and the ellipses close to a C point also vary out of the 
plane of circular polarization. As a way of understanding 
this, Freund [9] recently proposed polarization 'Mobius 
bands' around C points: on a loop enclosing the C point, 
the ellipse axis generically winds ±1/2 or ±3/2 times 
around the loop, this twist number being half-integer 
due to the fermionic C point index. Such Mobius band 



geometry is, however, very complicated, depending on 
both the in-plane and out-of-plane C point geometry and 
no mathematical analysis has been made of them. 

Here, I will describe an alternative, simpler geometric 
classification, the 'O classification' of the out-of-plane 
variation around a C point in 3D complex vector fields: 
the number of directions out from the C point in which 
the axis lies in the plane of the C point is generically one 
or three. The odd number is the signature of a fermionic 
multivaluedness: there are two opposite ellipse axis vec- 
tors, which must be exchanged in a circuit around the 
C point. This contrasts with the 'bosonic' out-of-plane 
variation around a typical polarization ellipse, for which 
there are generically two opposite directions in which the 
local ellipse axes lie in the the central polarization plane 
(as described below). Examples of the possible fermionic 
out-of-plane structures are represented in Fig. [1] 

In considering the mathematical basis of this O clas- 
sification, it is natural to choose a local helical basis 
E = (£7-|- , E_, E z ) for the complex vector field around 
the point of interest in three dimensions (placed at the 
origin), with E z = (i.e. the polarization ellipse is in 
the xy plane). In this basis, the helical unit vectors are 
e± = -4-(l,±i, 0) with respect to cartesian coordinates, 
and l-E+l > \E_ \ without loss of generality. Since E z = 
at the origin, there is a complex phase singularity in the 
local 2-component, and at an ordinary ellipse, E + ,E_ 
are nonzero constants. At a C point, E— = at the ori- 
gin as well, leaving only E+ nonzero, and so in this plane 
the C point appears right-circular. 

The choice of a helical basis makes the in- plane C point 
classification straightforward. The projected, in-plane el- 
lipse axis direction is given by half the argument of 
a = Si + 1S2 = E* + E- where Si, S2 arc two-dimensional 
Stokes parameters in the plane [TOllllj : at a C point, this 
in-plane axis angle singularity is directly related to the 
phase singularity in E_ . It is further convenient to choose 
local coordinates and phase such that Vi?_ = (u x , iv y , 0) 
with \v y \ > u x > 0, so the sign of v y determines the 
C point index, and the anisotropy ellipse ^2rtl4] of the 
phase singularity in E- (and thus the C point) is aligned 
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Fig. 1. (Color online) Representation of polarization axes 
close to C points, (a) E = (\/2 + x + iy, iy/2 — ix + 
y, (-1.40-.18i)a;+(-1.40+.18i)j/). (b) E = (1.-0.98x4- 
0.20h/, i4-0.98ix4-0.20y, 0.63(x-2h/)). In each frame, the 
black lines represent the ellipse axes around the C point 
(sphere), in its plane of polarization. The lighter (darker) 
gray ellipse is the anisotropy ellipse for £L (E z ), and the 
blue surface above represents ±a z around the C point, 
with either three zeros (a) or one zero (b), resembling a 
2-shcet Riemann surface. 

with the xy-axes. With this choice, — ^axgE + = (3 is 
the in-plane axis orientation angle along the local x- 
direction, and in Ref. |5] it was shown that the index and 
L type depends only on the angle /3 and the C point 
isotropy, parametrized by T = 2u x v y / (u* 4 v 2 ). 

The complex quantity E-E w E+E- close to a C point 
is related to a. \E ■ E\ = \a\ is a measure of the polariza- 
tion ellipticity, and \ arg E ■ E is the vibration phase, de- 
termining where KcE is around the ellipse |2l6lfT?j , which 
is zero when the real part is along the major axis of the 
ellipse, and undefined for circular polarization. With the 
coordinate and parameter choices above, at a C point, 
\ aigE ■ E = —f3 with respect to the +x-direction. 

The director field of ellipse axes is conveniently repre- 
sented by the smooth field via the vibration phase [TTIfTB] 

a = ±ReVE* • E*E. (1) 
At a C point, a = 0, and in a loop around the point, 
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Fig. 2. (Color online) Surface Do = for alignment 
9 = 0, in (T, T z , x)-space, where Do > on the red 
side. The surface is periodic in x> an d in the four faces 
where T,T 2 = ±1, the line Do = resembles the L 
classification locus [5]. 

± — > =p, revealing the singularity's fermionic nature. The 
O classification is determined by the z-componcnt a z in 
a neighborhood of the C point. 

Close to a typical point of elliptical polarization, 
placed at the origin with the coordinates chosen above, 

a z w ±Rc^E* + E*_(x,y,0) ■ VE Z , (2) 

where E+E*_ is a nonzero complex constant evaluated at 
the origin. Therefore a z = locally along the two lines 
perpendicular to ±Rey / E+El VE Z , and ±a z is repre- 
sented by a pair of tipped planes through the C point. 

At a C point, both E- and E z are 0, and in this case, 
the out-of-plane component of the axis, in the coordi- 
nates chosen above, is 

a z w ±Re\E + \ l ' 2 e 1 ^u x x - iv y y{x,y,0) ■ VE Z . (3) 

The variation of a z therefore depends on the shapes 
of the two phase singularities, one in _E_ (which is 
square rooted), and the other in E z . The structure may 
be understood by writing the complex vector VE Z = 
e 2w (cos 6Ux - i sin 0V Y , sin 0U X + i cos 9Vy ) , with \V Y \ > 
Ux > 0. This can be interpreted as the anisotropy el- 
lipse of the nodal point in E z having axes proportional 
to Ux , Vy (in the X, Y direction) , making an align- 
ment angle 9 with the x, y directions defined by the 
axes of the anisotropy ellipse of the node in E_. The 
isotropy parameter of this E z anisotropy ellipse is writ- 
ten T z = 2U X V Y /(U X + Vy), similar to T defined above 
for the in-plane isotropy. The angle v is the vibration 
phase associated with the complex vector VE Z : the 
phase factor multiplying the complex E z to find the real 
a z depends on a combination of v, the vibration phase 
of the in-plane pattern /?, and the alignment angle 9. 
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The O classification of C points, that is, whether 
a z = once or thrice on a small loop around the C 
point, is determined in a similar way to the L classifica- 
tion [3HI]- The loop (x, y) = 8{cost,sint\v y \/u x ) 1 cor- 
responding to a circuit on the anisotropy ellipse itself, is 
a more convenient choice than the circle. Setting a z = 
from Eq. ([3]) on this loop and rearranging gives a cubic 
polynomial equation in tan(i/2), with real coefficients 
depending on the in-planc isotropy T, the out-of-planc 
isotropy T z the anisotropy ellipse alignment angle 9, and 
the phase x = ft + a combination of the vibration 
phases of the two anisotropy ellipses. This polynomial 
has one or three roots, corresponding to one or three ze- 
ros of a z on the loop (01, 03 respectively), depending 
on the sign of the discriminant of this cubic, given by 

D = 2(3 - T 2 - UTT Z - T 2 + 12T 2 T 2 ) 

- (1 + T) v 7 ! -T 2 (l + T z f cos4( x - 0) 

+ 2(2 - T)(l + T)(l + T z )y/1 -T 2 cos2(2 x - 9) 

- 4\/l - T 2 (2 - 7TT 2 )v/l - T 2 cos 29 

+ 2(1 - T 2 )(l - T 2 ) cos4<9 - 6y/l - T 2 (l - T 2 ) cos4 X 
+ 2(1 - T)(2 + T)(l - T,)v/l-T2 C os2(2x + 9) 

- (1 - T) y/1 - T2(l - T 2 ) 2 cos4( x + 9), (4) 

with D > for 03 (< for 01). Equation (gj is the 
main result of this Letter. 

The out-of-plane classification is determined by Equa- 
tion (j4j), depending on the four parameters T,T Z ,6, 
and x^ which is therefore hard to visualize. A particu- 
lar slice through this four-dimensional parameter space, 
with 9 = (the two anisotropy ellipses are aligned), is 
shown in Fig. [2] The vertical faces of the figure demon- 
strates a property of the set Do = : when the in-plane 
C point is isotropic (T = ±1, corresponding to a uniform 
lemon or star) , or the out-of plane node in E z is isotropic 
(T z = ±1), then 9 is not defined and the Do = locus, 
depending on \ an( i T or T 2 , is a cusp, similar to the 
corresponding Dl discriminant determining the in-planc 
L classification, determined by T and (3 [8]. 

This new geometric classification, characterizing the 
out-of-planc behavior of a nonparaxial C point, com- 
plements the traditional index, line and other in-plane 
C point properties. Combined with the L classification, 
it also gives some insight into the Mobius index [9]. A 
polarization Mobius band may be thought of as a rib- 
bon around the C point with axis curve (5(cos <fi, sin <fi, 0), 
framed by the projection of a into the normal plane 
of the curve at each point. In order for this ribbon to 
have ±3/2 twists, this projection must lie in the plane 
a z = three times, and also be perpendicular to the 
curve (i.e. (a x ,a y ) tangent to the circle) three times; a 
necessary condition is that both Dl and Do > 0. 

The L and O classifications correspond to the count of 
special directions in the x, y plane around a C point. For 
the L classification, it is the number of times that (a z , a y ) 
is tangent to the circular loop (equivalent to the number 



of times it is radial). The actual positions of these lines, 
in fact, determine whether the Mobius index is ±3/2 
or ±1/2 precisely, depending on how these two sets of 
three directions interlace. If, in a circuit around the C 
point, the special directions alternate, then a must twist 
around the C point exactly 3/2 times, whereas, if not 
(i.e. there is a repeat of the same type of direction) the 
ribbon twist sense changes direction, resulting in ±1/2 
a twist. In principle, the two sets of three directions can 
be located precisely as roots of the corresponding cu- 
bic polynomials in tan <f>, although this method is alge- 
braically taxing, and further information is required as 
well to determine the sign of Mobius twisting. 

Beyond the O classification of C points, there is an 
equivalent out-of-planc classification for points of linear 
polarization [JUS], and their twist indices [T7], and it 
may be extended to account for more twists of larger 
loops beyond first order [18] - Furthermore, this new clas- 
sification is very general, although challenging to mea- 
sure in a real 3D optical field; similar kinds of topological 
defect, such as disclinations in liquid crystals, are likely 
to have similar variation (related to the Frank 'twist' 
parameter close to the defect Qjj]), exhibiting further 
geometric fermionic properties of topological defects in 
director fields. 

I am grateful to Isaac Freund, Robert King and Brina 
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